o 

(N 
3 



oo 
in 



O 
O 



Horizon Problem Remediation via Deformed Phase 

Space 

S.M. M. Rasouli* Mehrdad Farhoudi^ and Nima Khosravi* 
Department of Physics, Shahid Beheshti University, G.C., Evin, Tehran 19839, Iran 

April 26, 2011 



Abstract 



We investigate the effects of a special kind of dynamical deformation between the momenta 
of the scalar field of the Brans-Dicke theory and the scale factor of the FRW metric. This 
special choice of deformation includes linearly a deformation parameter. We trace the deformation 
footprints in the cosmological equations of motion when the BD coupling parameter goes to infinity. 
One class of the solutions gives a constant scale factor in the late time that confirms the previous 
result obtained via another approach in the literature. This effect can be interpreted as a quantum 



5-h , gravity footprint in the coarse grained explanation. The another class of the solutions removes the 

big bang singularity, and the accelerating expansion region has an infinite temporal range which 
overcomes the horizon problem. After this epoch, there is a graceful exiting by which the universe 
enters in the radiation dominated era. 
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1 Introduction 



■ i 

^ ■ General relativity (GR) and quantum theory, as generally believed, are two prominent paradigms 



to illustrate the nature. The idea of GR was inspired by the general covariance, the principle of 
equivalence and also the Mach ideas. However after the formulation of GR established by Einstein, 
its full satisfaction with the Mach principle has been a matter of debate. That is, although the 
matter content of the universe affects the geometry, but there are still vacuum solutions for GR, 
contrary to the strong version of Mach idea that states if there is no matter then, there will be no 
geometry. This proposition has made further considerations for alternative gravitational theories 
to be more Machian. One approach to this purpose has been stated by the scalar-tensor theories 
which among them, the Brans-Dicke (BD) theory is the simplest one [1]. In this scope, a scalar 
field plays the role of Newtonian gravitational constant and makes the theory to be more Machian^ 
Also in the BD theory, there exists an adjustable dimensionless parameter which, in principle, can 
be fixed by observations. On the other hand, with suitable boundary conditions, GR is deducible 
from the BD theory when making this parameter goes to infinity limit (though not always, see, e.g., 
Refs. [3, 4, 5, 6]). In this view, the BD theory presents a modified version of GR (at least for traceless 
energy-momentum tensor [6]). Hence in this sense, there may also exist some other features in the BD 
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1 However, the other motivation for development of this variant of GR was the Dirac large number hypothesis [2]. 
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theory in addition to those that are usually in GR. Such a kind of features (if any) can be interpreted 
as internal symmetries (or structures) of GR which are somehow missed if one just starts purely with 
the Ricci scalar. This aspect can also be viewed as if GR is only an effective theory, and the BD 
theory when the u goes to infinity is a more realistic one. 

The quantum theory is mainly established to justify the behavior of very small scales. Thus, it 
should be considered for small scale behaviors of GR (which is an excellent theory for large scale 
structures) as well. For example, in the standard cosmology, the universe has been commenced by 
a big bang in the very early universe. It is generally believed that, in the big bang, the scale of 
the universe is almost zero, thus it is predictable that the quantum behaviors to be significant in 
this regime. Indeed, there have been considerable attempts to combine quantum theory and GR 
in order to achieve a quantum gravity theory [7, 8, 9, 10]. One of these approaches is deformation 
of the phase space structure [11] that introduces, if not all but at least, a part of the quantum 
effects. This mechanism has been employed in the context of cosmology [12, 13], in affecting the 
small and even large scale behaviors, by removing the singularity and by the coarse grained effects, 
respectively. In this approach, usually a (length) parameter, which can be interpreted as the Planck 
(length) constant, presents the quantum regime. However, it is crucial to recover the standard results 
by taking the appropriate limits of this parameter. The another significant aspect of this kind of 
mechanism is its correspondence to the other methods of quantization, which is not only interesting 
but also indispensable, see, e.g., Refs. [14, 15]. 

As the scalar-tensor gravitational theories involve more degrees of freedom, they can give more 
number of solutions than GR [16]. Hence, we purpose to investigate deformation of the phase space 
structure in this context by studying the flat FRW metric and employing noncommutativity between 
the momenta of the BD scalar field and the scale factor. Then, we obtain the BD dynamical equations, 
however for simplicity and to be able to compare the outcomes with the corresponding results in GR, 
we solve them when the BD coupling parameter goes to infinity. Though, it is important to note 
that such a procedure does not make a precise transition to GR when the trace of energy-momentum 
vanishes [6], but we should also emphasize that if one just starts the procedure by the GR formalism, 
one cannot achieve such a wide classes of solutions, as expected. 

In the next section, we briefly describe the BD theory in the context of Hamiltonian formalism 
which is essential for introducing deformation (noncommutativity) in the phase space. In Section 3, 
we derive the BD equations of motion in the presence of a special kind of deformation, then deduce 
their cosmological implications when the BD coupling parameter goes to infinity. Then, we discuss 
the solutions for different signs of the deformation parameter and the integration constants of the 
solutions while highlighting their effects. Finally, we will end up the work by conclusions in the last 
section, while two appendixes have also been furnished. 

2 Hamiltonian Formalism for Brans— Dicke Theory 

In this section, we review the BD theory in the Hamiltonian formalism. However, in order to study 
the cosmological behavior, we consider the spatially flat FRW metric as the background geometry, 
namely 

ds 2 = -N 2 {t)dt 2 + a 2 {t) (dx 2 + dy 2 + dz 2 ) , (1) 

where N(t) is a lapse function and a(t) is the scale factor. 

The BD Lagrangian density in the Jordan frame [1, 17, 18] in vacuum is given by 

C = y/^^>R-^srM^j, (2) 

where the Greek indices run from zero to three. Also, the (j) is the BD scalar field, R is the Ricci 
scalar and uj is the BD coupling constant that is supposed to be bigger than —3/2 for non-ghost 



scalar field situations [19, 20, 21, 22]. Replacing the Ricci scalar associated to metric (1) into the 
above Lagrangian gives 

C = -6iV _1 ad 2 - 67V _1 a 2 d0 + ooN^a 3 ^ 1 <j> 2 , (3) 

where the dot represents derivative with respect to the time and a total time derivative term has 
been neglected. Thus, the corresponding Hamiltonian is 

no = ^ (-^a-V 1 ^ + ~ \a' 2 P a P,) , (4) 

where \ = 2uj + 3. As the momentum conjugate to N(t) vanishes, one has to add it as a constraint 
to the above Hamiltonian. Therefore, the Dirac Hamiltonian becomes 

n = n + xp N , (5) 

where A is a Lagrange multiplier. Here, we consider the ordinary phase space structure described by 
the usual ordinary Poisson brackets 

{a,P a } = {4>,P (/> } = {N,P N } = l, (6) 

where the other brackets vanish. Therefore, the equations of motion with respect to the Hamiltonian 
(5) are 

a = { a ,H} = -^(^a- l 4>- 1 Pa + \a- 2 P (p y (7) 
Pa = {Pa,n} = -^^a-^- 1 P 2 a -la'UPl + a~ 3 P a P^ (8) 
4> = {4,H} = -(a-^P (j> -l-a- 2 P a ) , (9) 



X V T " 2 

= {P h U} = -?L^a- 1 ct ) - 2 P 2 a+\a- 3 pi), (10) 

N = {N,H} = \, (11) 

Pn = {PN,n} = ^[^a- l ^ 1 P 2 a -] 1 a-^Pl+ l -a- 2 P a P^. (12) 

Let us work in the comoving gauge, that is we fix the gauge by N = 1. Also, to satisfy the 
constraint Pn = at all times, the secondary constraint Pjy = should also be satisfied. Hence, by 
Eq. (12), one obtains 

p ^ = K i± vl) arlPa - (i3) 

Now, differentiating Eq. (7) with respect to the time, while using Eqs. (8), (9) and (13), leads to 

?* ) a-'h 2 . (14) 

In addition to the trivial static solution, one can get solutions as a(t) = C\(t — t- m i) q± where C\ and 
the initial time ti n j. are integration constants, and q± is 



Q± 



3x-l 



(15) 



when x 7^ 1/3 (i.e. u / —4/3). In what follows, we assume that t; n i. = 0, and the constant C\ can 
be fixed by the scale of the universe at an appropriate definite time. Also, one can easily obtain 

j = T f zyzx \ a ; (16) 



with solutions <fi(t) = C2t s± , where C2 is an integration constant and s± is given by 

when again \ 7^ 1/3- Through the Hamiltonian approach, we have actually rederived the O'Hanlon 
and Tupper solution [23, 24], as expected. This solution has a big bang singularity when t tends to 
zero. Also, we should emphasize that the behavior of the scale factor and the scalar field depends on 
the BD coupling parameter. In particular, when co goes to infinity, the non-trivial solution 2 is not 
the same as the corresponding result of GR (i.e. the Minkowski space-time), for in this case one gets 

a(t) = Cit 1/3 and 4> = constant, (18) 

where the scale factor has a decelerated expanding behavior. 

In the case x = 1/3, the solutions for upper sign are a = C^(t — t' ) 2 ^ 3 and <f> = C±{t — ig) -1 , 
and for lower sign are a = aoexp^/ig) and 4> = 4>o exp(— 3t/£o), where C3, C4, do, (fro, t' and are 
constants. 

In the next section, we investigate a modified version of the above formalism by introducing a 
noncommutative model. 



3 Deformed Phase Space Brans— Dicke Structure 

As mentioned in the introduction, deformation of the phase space can present a sort of tracing 
quantum footprints in a given model [11]. Of course, a general deformation makes equations very 
difficult and even unsolvable. Hence, it is customary to pick a proper choice which not only makes 
calculations be possible, but also gives non-trivial results. Indeed, by appealing to the simplicity 
principle (or the Occam's razor), simplifications are usually performed in most toy models, and even 
in real ones, before a consistent and complete theory is deduced. Thus in this work, we just consider 
a dynamical deformation between the conjugate momentum sector as 

{P' a ,P' (t> } = l<t>'{t) (19) 

and leave the other Poisson brackets among the primed parameters (corresponding to those appeared 
in relation (6)) unchanged. Hence, the Jacobi identity is still satisfied. In general, noncommutativity 
between the momenta is a kind of generalization of the usual noncommutativity between the spatial 
coordinates [25]. In addition, noncommutativity between the momenta, in effective, has similarity 
with the behavior of a charged particle in the presence of a magnetic field. In the scope of gravitational 
theories, this kind of noncommutativity can be interpreted as a gravitomagnetic field [26]. Also, one 
may find a clue in the string theory, especially in the flux compactification, as mentioned in Ref. [27]. 
The dynamical behavior of deformation as time dependence has also been employed in the literature, 
e.g. the K-Minkowskian spacetime [28, 29] and the generalized uncertainty principle [30, 31], that 
are considered also in the cosmological phase space [32, 33]. In addition, the main reason for the 
peculiarity of the chosen term in the right hand side of (19) is the dimensionality analysis that is 
described in more details in the Appendix A. 

The minimally deformed (noncommutative) version of Hamiltonian (4) is achieved by replacing 
the unprimed variables with the primed ones, namely 

K = ^ (-^a'-V'-^ 2 + VV^ 2 - \a'- 2 P' a P',) . (20) 

However, it is more convenient to re-introduce the new variables by applying the standard transfor- 
mation [25, 34] 

1% = P^ - lac/>, (21) 

2 One could easily obtain the solution from Eq. (14) for this limit, however for the sake of completeness, we have 
firstly derived explicit solutions for any oj too. 



A 



where the other unprimed variables are equivalent to their primed counterparts. By considering 
the above transformation, relation (19) is satisfied when the unprimed variables satisfy the ordinary 
Poisson brackets. Then, substituting P^ from (21) into Hamiltonian (20) gives 

?C = 7io + — " a"V + , (22) 

where we have also substituted H' , as a function of the primed variables, with Hq c , as a function 
of the unprimed ones via the employed transformation. Once again, the noncommutative Dirac 
Hamiltonian is 

U ™ = U ™ + \P N . (23) 
Therefore, the equations of motion become 

a = {a J U ac } = ^^a- l r 1 Pa-\a- 2 P 4 ,+ l -a~ 1 ct>j, (24) 

Pa = {Pa,n nC } = 

~ [^^-'Pa ~ \a- A <t>Pl + *- 3 PaP<f> ~ £a"V + 2la~^P, - l -a-^ , (25) 
= {<P,n nc } = ^{a- 3 ^P ( f > -^a- 2 P a -la- 2 ^y (26) 

P, = {P 4 ,,Hn = ~(^a- 1 r 2 P a 2 + la- s P 2 + ll 2 a-^ (27) 

N = {N,H nc } = \, (28) 
P N = {P N ,H nc } = 

+1 feci- V" 1 ^ " ^a- 3 ^ + " ^a" V 3 + IcT^P* - ^P a ) . (29) 



In the comoving gauge, i.e. N = 1, the secondary constraint Pat = gives 

a 



P *=2 



2l4> + (1 ± 



Employing Eqs. (24)-(27) and (30), and performing a little manipulation lead to 



-\ \ -1.2 , 1 ,..-2 



and again 



(30) 



a = — ;= a a ± . /a <*d (31) 



4 = T f * (32) 



f=T-|^=. (34) 



Solution of Eq. (32) can be in the form 

4> = <j)o a 5 , (33) 

where 0o is a constant and £ is 

Substituting (33) into (31) yields 

d = -( 2 * ) a^d 2 ± /a«- 2 d. (35) 

Note that, when the deformation parameter I tends to zero, all noncommutative equations reduce to 
their corresponding ones in the previous section. 



Now, as proposed, we are interested to investigate effects when uj goes to infinity, though again it 
does not mean that it makes transition to the standard GR, as has been shown for the commutative 
case in the previous section. However, as it is obvious from Eq. (35), in order to obtain such effects, 
it crucially depends on how the constant <po is, or can be, related to the uj. Actually for this purpose, 
the value of constant </>o (which represents different initial conditions) proportional to yfx can be a 
reasonable one. A particular motivation for it, however, is the new term in Eq. (35) (in comparison 
to Eq. (14)) which should not vanish in the limit \ ~ * 00 on t ne one hand. Namely, if it would 
vanish, one would not be able to see any new effects in comparison to the original undeformed theory 
in the limit uj oo. On the other hand, the new term also should not become infinite. However, a 
consequence of this choice is that when uj tends to infinity, then 0o goes to infinity as well. Though, 
this brings a technical problem, i.e. it makes some ambiguities in the behavior of 0— field in solution 
(33) when 0o — > oo and £ — > 0. In this limit, <p is a time independent (constant) variable which is 
infinity. To study the inconvenience caused by this divergence of 4>o, the renormalization argument 
may assist in the following manner. 

It is well-known that the procedure of renormalization occurs in the quantum field theoretical 
level. However in our toy model, this may indicate itself naively only in the deformation parameter 
as the only presenter of quantum regime in this work. Hence, as the first option, the deformation 
parameter I = /bare can be re-defined in an appropriate way that makes the transition from Eq. (35) 
to Eq. (36) being possible. That is, it can be re-defined as Prenormalized = 4>o ^bare/V^X ( see below 
Eq. (36)) with a finite 4>o. Then, when x — > 00 , the l\, are deformation parameter goes to infinity such 
that the Prenormalized deformation parameter becomes a finite constant. Also, there is an alternative 
approach which is considered in the Appendix B. 

Therefore, by taking uj goes to infinity and choosing the minus sign 3 in Eq. (35), one gets 

o = -2a~ x a 2 - £aT 2 a, (36) 

where <fiol = yJ12\£, which fixes the new parameter £ with dimensionality L~ l . Substituting a 2 a = 
(a 2 a) — 2aa 2 into Eq. (36) gives 

(a 2 a) = -ia, (37) 

that yields a 2 a = —£a + C, where C is an integration constant with dimensionality [C] = L _1 . Then, 
one easily obtains 



1 2 C C 2 C 
-a + —a + — 7T hi 



2 . p — , «<-« + *.), (38) 

where to is an integration constant too. Obviously, the above equation is invariant under the transfor- 
mation (£, C, t) — > (—£, —C, —t). This symmetry makes a counterpart relevant between the solutions 
and, consequently reduces the number of investigations for different cases by half? Thus, in the follow- 
ing categorization, we consider the two probable options (the Case I and Case II) of the logarithmic 
term in Eq. (38) only for interesting cases of different signs of the I and C, without probing the 
counterpart solutions. Also, for the sake of completeness, we explicitly investigate the solutions when 
i tends to zero in Case III. 

3.1 Case I: a - f > 

As mentioned, we investigate this case for different signs of the £ and C. 

3.1.1 Case la: Negative £ & C 

For convenience, assume £ = — £ > and b = —C > 0, thus Eq. (38) reads 

^a 2 + tfa + 2 In (a - 0) = £(t -to), (39) 

3 This choice is not restrictive, for in the following we will consider different signs for the I. 

4 For example, the case I < and C < is the counterpart of the another case £ > and C > when t — > —t. 



a 



where $ = C/£ = b/£>0 with valid domain 5 a > This means that the initial value of the 
scale factor cannot be zero and indeed, the universe has been started with a non-vanishing size. 
Therefore, one may interpret that the existence of a deformation parameter, as an indicator which 
usually presents the quantum corrections to models, removes the big bang singularity. Actually, this 
result is a common expectation in the quantum cosmological models. 
Then, by differentiating (39), for a > we obviously get 

a = -^(a-tf)>0 (40) 

and hence, 

I 2 

a = -_( a -0)(a-20). (41) 

Thus, the sign of a depends on two different regions, d < a < 2$ and a > 2$, which we investigate in 
the following. 

Region -d < a < 2-d 

In this region, the expansion is an accelerated one in similar to the inflationary phase. Though, 
it is not exactly as the standard inflationary phase, but it can solve the horizon problem as will be 
discussed in the following. Usually, a successful candidate for the standard inflation should satisfy two 
essential properties among the other ones, namely the 60 e-fold duration and a graceful exit from this 
epoch. Our model naturally satisfies the latter requirement, for the scale factor transits to the next 
region, i.e. a > 2-d, where it decelerates. However, at first glance, it looks that it does not satisfy the 
former condition as it has much less than 60 e-fold duration. Indeed, the number of e-fold definition, 
i.e. N = afinai/ainitiai, for our model is N = 2-d/d = 2. Nevertheless, its result is comparable to the 
standard inflation one by presenting a solution for the horizon problem, which we indicate it after a 
brief review on the successes of the standard inflation while clarifying the horizon problem. 

It is well-known that the most important problem of the standard cosmology, which is solved by 
proposition of an inflationary scenario, is the horizon problem. Of course, the inflation also solves 
the relic particle abundances (or the monopoles) and the flatness problems. However, it is generally 
believed that among these problems, the horizon problem is the most important one, for at least 
there are alternative scenarios that can resolve the other two problems in the same manner as the 
inflation does [35, 36]. The horizon problem arises when the universe is observed to be isotropic 
and homogeneous in the large scale structure. This requires that the initial conditions must be in 
a way which give such a universe. The problem with the standard cosmology is that although the 
matter fluctuations have been inhomogeneous at the initial level, but these fluctuations did not have 
enough time for interactions and transforming information about their situations. Consequently, 
the inhomogeneous initial conditions should naturally result in inhomogeneous present large scale 
structure which is in contradiction to the observations. The inflationary idea solves this problem 
by taking a homogeneous part of the initial condition, and inflates it to an appropriate size for the 
beginning of the radiation dominated era. In our model, the horizon problem is solved in another 
way. 

Actually, the accelerating phase occurs during a; n iti a i — > "& and afi na i = 2i9, which is from 
^initial — > — °o to a finite final time, that is tfi na i = "9 2 (4 + In d)/£ (assuming to = 0). This means that 
the accelerating phase takes infinite time, At = (ifi na i — initial) — ► °°> an d during this phase, the 
matter fluctuations can interact with the other parts of initial conditions, exchange information about 
their local structures, and hence, approach to an equilibrium state which is presented by a homoge- 
neous structure. Thus, although in our simple model, the accelerating phase cannot be interpreted 
as the standard inflationary era, but it can address the horizon problem of the standard cosmology. 

5 We have neglected the equality a = for it makes t becomes — oo. 
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Fig. 1: The solid line shows the behavior of the scale factor for Case la with i — C = —0.1. Below the dashed line 
(i.e. ■& < a < 2$ — 2C/£ = 2), one has an accelerating phase, and above it (i.e. a > 2i9 = 2C/1 = 2), a decelerating 
phase. In the right figure, the dotted curve represents the a(t) which is negative for a > 2i9 = 2. Note that, the a(t) 

curve has been rescaled for a better clarification. 



This quasi static 6 accelerating phase is very similar to the Hagedorn phase of string gas cosmology [37]. 
Region a > 2-d 

In this case, the expansion of universe is decelerating, and when a — > oo, the first term in 
Eq. (39) is the dominant one, hence in this limit, a(t) tends to i 1 / 2 that behaves as the radiation era. 
This phase occurs exactly after the above accelerating phase, and can be interpreted as the radiation 
dominated phase after the usual inflationary epoch in the standard cosmology. Indeed, this result 
is completely in agreement with what is usually proposed for the universe in different cosmological 
models, the standard cosmology with or without an inflation. 

The behavior of the scale factor for Case la is plotted as a solid line in Fig. 1. 

3.1.2 Case lb: Negative I & Positive C 

This case is very similar to the previous one, except the sign change in the acceleration. Hence, 
for the entire valid region of the scale factor, i.e. a > 0, the a > and a < give a decelerated 
expanding behavior. When t — > oo then a(t) tends to i 1//2 , exactly as the previous case, however here, 
there are not interesting properties. The behavior of the scale factor is depicted in Fig. 2(left). 

3.2 Case II: a - f < 

Once again, we investigate this case for different signs of the I and C too. 

3.2.1 Case Ha: Positive I & C 

To describe this case, let us repeat Eq. (38) as 

1 



^a 2 + i?a + tf 2 In (0 - a) = i(-t + t ), 



(42) 



where the valid domain of the scale factor is < a < 0, which the a > is a physical constraint. 
That is, in this case, there is no regulator to prevent the big bang singularity, for a = occurs at 



3 For it takes infinite time to double the value of the initial scale factor. 
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Fig. 2: The solid line in the left figure shows the behavior of the scale factor for Case lb with £ = — C = —0.1. As it is 
obvious, the scale factor asymptotically is proportional to t 1 ^ 2 (the dotted line). The solid line in the right figure shows 
the behavior of the scale factor for Case Ha with £ = C = 0.1. When t — > oo, the scale factor approaches ■& — C/£ = 1. 

t = t - ($ 2 ln$)/t. However, the upper bound of the scale factor approaches to a constant when t 
goes to infinity, i.e. for the late time, one gets 



This approaching to a constant value is a direct consequence of the existence of a deformation pa- 
rameter, and more interestingly, it shows itself far from the big bang. As mentioned before, since this 
deformation parameter may be interpreted as a consequence of the quantum effects, then this feature 
may also be viewed as a quantum gravity effect when the scale of the universe is significant. That is, 
this behavior can be a phenomenological property for the quantum gravity. However, these kind of 
modifications can be perceived as a semi-classical model or, as a model beyond the BD theory but 
still in the classical regime. The scale factor decelerates in this choice, and its diagram is plotted in 
Fig. 2(right). The graph shows that the behavior of the scale factor is in agreement with the results 
obtained by Ref. [32] in where a dynamical deformation between the lapse function and the scale 
factor has been employed. 

3.2.2 Case lib: Positive £ & Negative C 

The case I > and C < gives negative values for the scale factor which is not acceptable. 

3.3 Case III: £ Tends to Zero 

As mentioned before, all equations reduce to their corresponding commutative ones when the 
deformation parameter vanishes. Now, let us explicitly investigate it for Eq. (38). Hence, taking the 
limit I — > for the logarithmic term in Eq. (38) gives 



where higher order terms in I can be neglected. And obviously, the second and the third terms in 
the above relation cancel the second and the first terms in Eq. (38), respectively. The first term in 
relation (44) is a constant and can be absorbed by re-definition of to- Consequently, the scale factor 
tends to i 1//3 which recovers the commutative solution (18) when u goes to infinity, as expected. 




(43) 



C 2 ,C \ C 2 C Ca a 2 (a 3 



(44) 



n 



4 Conclusions 



We have introduced a deformation in the phase space structure of the two existing fields of the 
BD theory in the spatially flat FRW metric. Also, we have traced the quantum footprints in the 
cosmological equations of motion in the comoving gauge. All the noncommutative equations are 
shown that do reduce to their corresponding counterparts when the deformation parameter tends 
to zero, as expected. Then as proposed, we have investigated the effects when the BD coupling 
parameter goes to infinity. In this process, we have faced an integration constant that depends on the 
initial conditions, however in order to be able to trace the effects, we assume it to be proportional to 
the square root of the BD coupling parameter. We have also discussed our justifications for why we 
have to fix it in this way and to render other side effects. 

Finally, different cosmological results have been deduced due to the different possible signs for 
the two arbitrary parameters of the solutions, namely the deformation parameter, I, and the another 
integration constant (i.e. the C in Eq. (38)). For one class of the solutions, the result predicts a 
constant value for the scale factor in the late time that is in agreement with the results obtained in 
Ref. [32]. This feature may be interpreted as a quantum gravity footprint in the large scale. 

A more interesting result is achieved by the another class of the solutions. In this case, it is shown 
that the existence of the deformation parameter (or equivalently, the quantum correction) removes 
the big bang singularity by preventing the scale factor tends to zero. This also has an infinite temporal 
range for an accelerating expansion region. However, this phase is not the standard inflationary phase, 
for its e-fold duration is a very small number, but it can appropriately overcome the horizon problem 
that is the main one in the standard cosmology. Indeed, due to the very long time duration of this 
phase, the matter fluctuations can transmit their information to the other parts of the universe, and 
consequently become homogeneous. Implicitly, after this epoch, there is a graceful exiting and then, 
the universe enters in a radiation dominated era which is naturally in agreement with the standard 
cosmology, with or without an inflation. Note that, these consequences are held just by introducing 
a constant parameter without considering any potential in the model. It should also be mentioned 
that the model just makes a (classical) background plausible to address the horizon problem similar 
to (classical) background of the inflationary models. 

However, one of the major success of the standard inflation is its prediction of (quantum) fluc- 
tuations' behavior. The standard inflation anticipates a scale invariant spectral index which is in 
good agreement with the observations. For our model, considering the fluctuation of dynamics is 
important, for not only to compare with the observational data but also, to test the stability of the 
model. It should be checked whether inhomogeneous arbitrary initial conditions can have a significant 
effect in the late time behavior or not. To overcome such general questions, one needs to perform 
more investigations, perhaps employing the perturbative analysis which in our model is still more 
complicated than in the standard inflation, due to the existence of the BD scalar field as well as the 
deformation parameter. These are interesting investigations for further considerations, and are not 
in the scope of the current work. 

In fact the above achievement in solving the horizon problem can be viewed as a natural con- 
sequence of the noncommutativity approaches. That is, it is well-known that the varying speed of 
light and noncommutative models are related to each other [38], where the first motivation for the 
former models has been raised in order to achieve an alternative approach to the standard inflation. 
Indeed, the coordinate noncommutativity has been employed in the cosmological context for the same 
purpose as well, see, e.g., Refs. [39, 40]. 

Appendix A: On Dynamical Deformation As Relation (19) 

Let us first indicate the dimension of the Poisson bracket {P a , P(j>}- In this work, we have employed 
the units H = 1 = c, therefore, from the Plank length, lp = y/fiG/c 3 , the dimensions of G and <fi 
are [G] = L 2 and [4>] = L~ 2 . The scale factor and the lapse function are dimensionless parameters, 



1 n 



the dimensions of coordinates and the BD Lagrangian are [x^] = L and [C] = L~ 4 . Hence, one can 
conclude that [P a ] = L~ 3 and [P^] = L^ 1 , and consequently [{P a ,i-^}] = L^ 1 . On the other hand, 
the dimension of the deformation parameter is [I] = L. Therefore, from dimensionality aspects of 
view, the dynamical deformation (19) is a plausible choice. Besides, from simplicity point of view, 
with this choice, no other extra field has been introduced in the model. 

Of course, one may also propose other choices that still can satisfy the dimensionality of {P a , P<f>}, 
but the suggested relation (19) is a first order (linear) term in the deformation parameter as well. 
This suggestion is also a length indicator that can present and trace the quantum behaviors, and if 
the length indicator vanishes, one will recover the standard (classical) counterpart relations. On the 
other hand, in cosmological models a length parameter (e.g. the Planck length that is a function of 
h) is physically a more realistic and plausible choice. That is, among different choices that can be 
selected as a quantum indicator, a length scale is an appropriate one for cosmological models, that 
can compare different scales for the quantum or classical aspects as well. 



Appendix B: Discussion on Fixing Integration Constant <pQ 



To fix the inconvenience behavior of the 0o when x ~ > oo, one may also apply the renormalization 
procedure alternatively to the matter field ^matter- However in our model, there is no matter field and 
G (which is equal to 1/0 in the BD theory) does not appear in the equations of motion and makes 
its value non-effective, but the above procedure can be applied as if the matter field is turned on. 
Hence, in this case, the value of Newtonian gravitational constant can be regularized by re-definition 
of the matter field as well as the BD scalar field. To be more specific, in the presence of a matter 
field, the BD Lagrangian is 



-fj 



<f>R- V>,M^ + ^matter) 



(B.l) 



Now, as an overall constant has no role in the form of equations of motion, one can multiply the 
above BD Lagrangian by a dimensionless parameter (</>oGWe) 1 to get 



/ 7, ( AT} ^ ± ± I t bare | 



00 Chare 



R-u 



00 Chare 







0oGbarc ' ,(i ^ ^oGbarc / ,v 0oGWe 



+ 



-L 



bare 
matter 



(B.2) 



The 00) that goes to infinity when oj — > oo, can be absorbed in re-definition of and ^matter ^ a 
renormalization process such that 

^•renormalized r 



UJ 



renormalized 



0i? - -r<r 0, M 0,, + l™ 



(B.3) 



where = G' r ^ ormalized = ^normalized = 0/(0 O G barc ) and L-Xr lized = ^e er /(0 oGbare ). Thus, 
the difficulty of infinite value of 0o can be solved, at least naively, by the renormalization procedure. 
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